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Abstract 

We study the Higgs sector of the U(l)'-extended MSSM with CP violation. This is an extension 
of the MSSM Higgs sector by one singlet field, introduced to generate the fi term dynamically. We 
are particularly interested in non-standard decays of Higgs particles, especially of the lightest one, 
in the presence of CP violating phases for fj, e ff and the soft parameters. We present analytical 
expressions for neutral and charged Higgs bosons masses at tree and one-loop levels, including 
contributions from top and bottom scalar quark sectors. We then study the production and decay 
channels of the neutral Higgs for a set of benchmark points consistent with low energy data and 
relic density constraints. Numerical simulations show that a Higgs boson lighter than 2mw can 
decay in a quite distinctive manner, including invisible modes into two neutralinos (h — > X°X°) U P 
to ~ 50% of the time, when kinematically allowed. The branching ratio into h — > bb, the dominant 
decay in the SM, is reduced in some U(l)' models and enhanced in others, while the branching 
ratios for the decays h — > t + t~, h — > WW* and h — > ZZ* — > 4£ are always reduced with respect 
to their SM expectations. This possibility has important implications for testing the U(l)' model 
both at the LHC and later at the ILC. 

PACS numbers: 12.60.Cn,12.60.Jv,14.80.Ly 
Keywords: Supersymmetry, Higgs, LHC 
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I. INTRODUCTION AND MOTIVATION 



Confirmation of the Higgs mechanism of the Standard Model (SM) of particle prysics 
demands discovery of the elusive Higgs boson, likely seen at ATLAS and CMS 2j at a 
mass around 126 GeV. The Minimal Supersymmetric extension of the SM (MSSM), which is 
arguably the best motivated extension of the SM, offers stabilization of the Higgs mass, and 
moreover agrees well with the SM predictions in certain portions of its restricted parameter 
space. For instance, for the upper limit of nih ~135 GeV of the MSSM h — > bb is the 
dominant decay mechanism (~ 60%) in the SM and in the MSSM. On the other hand, 
in gauge and Higgs extended supersymmetric models, the properties of the Higgs bosons 
can be substantially different from that of the standard supersymmetric model predictions. 
For instance, the addition of one singlet field to the MSSM Higgs sector provides new tree- 
level contributions to the F- and D-terms, which stabilize the Higgs mass naturally at a 
larger value [3| . While many models predict a light Higgs boson around the weak scale (say 
~100 GeV), it will take some time to differentiate whether the boson discovered at the LHC 
belongs to the SM gauge symmetry, its minimal supersymmetric version (MSSM), or even 
to another extension such as the gauge extended versions of the MSSM. 

Extensions of the gauge symmetry by an extra U(l) factor (supersymmetric or not) are 
arguably the simplest extensions of the minimal model. The best justification for these 
extended models arises from assuming grand unified theories of strong and electromagnetic 
interactions (GUTs). In GUT symmetries, it seems difficult to break most scenarios directly 
to SU(2)l x U(1)y, as most models such as SU(5), 5*0(10), or E e involve an additional 
£7(1) group in the breaking. In supersymmetric U(l)' models [4] (referred to as U(l)' models 
from now on), the number of the neutral Higgs bosons is increased by an additional singlet 
field (S) over that of the MSSM, and the vacuum expectation value (VEV) of the singlet (S) 
is responsible for the generation of the u term, which allows Higgs fields to couple to each 

n n, 

other |5|, |6|; while number of charged Higgs bosons in the U(l)' extended models remains 
the same as in the MSSM. The interest in the Higgs sector of the U(l)' models also comes 
from the fact that such models arise naturally from string inspired models jf], 7-[s|, or as the 
dynamical solution to the [i problem in gauge-mediated supersymmetry breaking (GMSB) 



10]. While in the MSSM and in the U(l)' models lightest neutralino is the best candidate 



for a LSP, for the latter the LSP is less constrained. 



In these models, the lightest Higgs boson could potentially behave differently from the 
SM or the MSSM Higgs boson due to its singlet nature. While a Higgs boson of mass rrih ~ 
126 GeV can be predicted by the SM, or by the MSSM, or by numerous other models, the 
coupling of the Higgs to the known fermions or bosons is not the same in all these models. 
This fact can be extrapolated not only from the number of the Higgs bosons but also from 
their production and decay mechanisms. 

Of all the Higgs bosons in a model, the properties of the lightest neutral state are the most 
interesting, also given its likely discovery already at the LHC An interesting possibility is 
that its decay could be partially into invisible modes (a possibility hinted at by the reduced 
branching ratios into fermions at the LHC), or that there is another Higgs boson lighter 

n n 

than the one at 126 GeV, which decays completely or almost so, invisibly |11| . |12| . This 
scenario is motivated by global fits to the data at the LHC which indicate that a Higgs 
boson branching ratio of 64% is still unaccounted for [13|. 

In SM the Higgs can decay invisibly only into neutrinos, and this branching ratio is 
< 0.1% [14]. A light Higgs boson with substantial branching ratio into invisible channels 
can occur in a variety of models including scenarios with light neutralinos, spontaneously 
broken lepton number, radiatively generated neutrino masses, additional singlet scalar(s) 
and/or right handed neutrinos in the extra dimensions of TeV scale gravity. Among these 
possibilities, invisible decay of the lightest Higgs into light neutralinos is interesting since 
the light neutralinos are well motivated candidates for the Lightest Supersymmetric Particle 
(LSP), providing viable relic density explanations^. Decays into light neutralinos are possible 
in models with non-universal couplings, where LEP limits can be circumvented [17J, and in 
models with a light dark matter candidate. For instance, a study [18] indicates that this 
is a possibility in Eq, where the lightest Higgs boson of the Exceptional Supersymmetric 
Standard Model E 6 SSM can decay into the lightest neutralino pairs more than 95% of the 
time (l9fl. 

Additionally, the Higgs sector in extended models could provide potential sources of CP 
violation beyond the phase of the CKM matrix, also important for the observed baryon 
asymmetry of the universe. These phases can affect the masses and couplings of the Higgs 
bosons to the gauge and matter fields of the model, as was shown in studies of Higgs sectors of 



the MSSM 20] and next-to minimal supersymmetric models (NMSSM) 2l|. The phases can 
* Note that in U(l)' models the LSP can be the singlino 
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also affect production and decay rates patterns, as we will show in this study. In this work, 
we analyze the mass spectra of all the Higgs bosons, and the production and decay rates 
(visible and invisible) for the lightest Higgs in the U(l)' extended form of the MSSM with 
CP- violating phases. The masses of Higgs bosons in the U(l)' with CP violating phases 
model have received attention previously [22j, but we include them here, for consistency 
with the determination of their decay properties. Thus, we re-visit the Higgs sector of U(l)' 
models and calculate the masses, and in doing this, we improve on the previous calculation 
by including contributions from both (s)top and (s)bottom sectors at one-loop level, and 
add the constraint that the lightest neutral state should have mass ~ 125 GeV. 



Motivated by the above considerations, we study anomaly-free U(l)' models to probe 
their peculiar Higgs sector consistent with the known (astrophysical and collider) bounds, 
which are included in our benchmark points. We add the scalar quarks and neutralino con- 
tributions, and calculate a complete spectrum for the latter, and insure agreement with the 
relic density, assuming that the lightest neutralino is the LPS . We then study the produc- 
tion and decay modes of the lightest neutral Higgs boson, with the purpose of unraveling 
the existence and consequences of invisibly decaying Higgs bosons within the U(l)' model. 



The outline of our study is as follows. In the following section (Section [IT]) we introduce 
our effective U(l)' model, with particular emphasis on the Higgs sector. We present tree- 
level ( III Ap and one loop mass evaluations ( III Bp . and then an analytical calculation of the 
charged and neutral Higgs masses (III Cp . We then introduce the neutralino spectrum (III DP of 
the U(l)' model, which contains two additional neutralinos from the MSSM. We include the 
constraints on the particle spectrum coming from low-energy measurements of CP violation 
in Section [Till in particular from electric dipole moments (IIII Ap and Ek (1III Bj) . Following 



the exposition of the model and its constraints, we present our numerical investigations 
in Section IIVI in particular for the lightest neutral Higgs boson production and decay in 
IIV Al comment on the second lightest neutral state in IIVBI We summarize our findings and 
conclude in Section |Vj The full form of analytical solutions for the masses can be found in 
the Appendices. 



II. THE U(l)' MODEL WITH CP VIOLATION 



We review here briefly the U(l)' model, with particular emphasis to the Higgs and the 
neutralino sectors, as these are relevant to our study. The superpotential for the effective 
U(l)' model is 

W = Y S SH U ■ H d + Y t U c Q ■ H u + Y b D c Q ■ H d , (1) 

where we assumed that all Yukawa couplings except for Y t and Y b are negligible. As can be 
seen from (pQ), by replacing the fi parameter with a singlet scalar (S) and a Yukawa coupling 
(Yg), we resolved the fi problem of the MSSM [6j; fi is generated dynamically through the 
VEV of the S field (see III Ah and is expected to be of order of the weak scale. 

In addition to the superpotential, the Lagrangian includes soft supersymmetry breaking 
terms containing additional terms with respect to the MSSM, coming from gaugino masses 
M a (a = 1, 1', 2, 3) and trilinear couplings As, A t and A b as given below 

- C soft = (J2 M ^aK + A S Y S SH U ■ H d + A t Y t U c Q ■ H u + A b Y b D c Q ■ H d + h.c.) (2) 

a 

+ m 2 u \H u \ 2 + m 2 d \H d \ 2 + m 2 s \S\ 2 + M 2 ~\Q\ 2 + M 2 ~\IJ\ 2 + M 2 ~\D\ 2 + M\\E\ 2 + Mf |Z| 2 . 

Using Renormalization Group Equations (RGEs) these soft SUSY breaking parameters are 
generically non-universal at low energies. However, in our numerical studies, we choose 
not to deal with the evolution of the RGEs and instead assign them values which do not 
contradict with the current collider bounds. As we are interested in CP violation, we assume 
some of the soft breaking terms to be complex, selected as the trilinear terms (A t fi,s) and 
the VEV of the Higgs field S, as these assignments do not conflict with present low energy 
data. 



A. The Higgs Sector at Tree-level 

The effective U(l)' model inherits two Higgs doublets H u , H d from the MSSM, and has 
an additional singlet field S, all of which can be expanded around their VEVs as 



\ e i8d / v d + (f) d + iip d 



(Hu) = ~ 7n \ u , (H d ) = —= . 

V 2 \v u + (f) u + iip u J V2 \ s/2H d 



e i9. 



(S) = —(v s + ci>s + vps) , (3) 



in which v 2 = v 2 + v\ = (246 GeV) 2 . The fields in the superpotential are charged under the 
U(l)' gauge group with charges Q, required by gauge invariance to satisfy: 

Qh u + Q Hd + Qs = } Q Q[i + Q U3 + Q Hu = 0, Q Q3 + Qd s + Qn d = 0. 

The effective // parameter is generated by the singlet VEV (S), defined as 

fjieff = [ie %e % where [i = -^f, (4) 

so that with this convention \x is always real. For the remaining parameters we adopt the 
convention that the parameters are real, and explicitly attach CP violating phases where 
needed. Explicitly, arg(Af) = 9 t and similarly 8^ refers to the argument of A\>. In order to 
differentiate the phase of As from that of S we use small and capital letters: arg(5*) = 9 S , 
axg(As) = 6s- For the Higgs fields, we assume 9 U = 9^ = to avoid spontaneous CP 



breaking (SCPV) in the potential, associated with a real CKM matrix [23], which conflicts 
with experimental observations. However, to keep our considerations as general as possible, 
one can also define a new phase 

# s = arg (H u ) + arg (H d ) + arg (S) =9 u + 9 d + 9 s . (5) 



A detailed analysis of the Higgs sector with CP violating phases is available in 22] and 
references therein, but it is sufficient to mention that we assume 9 S ^ 0, which in a more 
general context could be replaced by #s 7^ 0. The tree level Higgs potential of the effective 
U(l)' model is a sum of F-terms, D-terms, and soft supersymmetry breaking terms: 

Vtree = V D + V F + V soft , (6) 

where the terms Vd, Vf and V S oft are: 

Vb = ^{\H U \ 2 - \H d \ 2 f + 9 4(\H u m^ ~ \H U ■ H£) + 9 Mq u \H u \ 2 + Q d \H d f + Q S |S| 2 ) 2 , 

o Z Z 

V F = \Y S \ 2 [\H U ■ H d \ 2 + \S\ 2 (\H U \ 2 + \H d \ 2 )] , 

Vsoft = m 2 Hu \H u \ 2 + m 2 Hd \H d \ 2 + m||^| 2 + (A S Y S SH U ■ H d + h.c), (7) 

where the coupling constant g 2 = g\ + g\- For the numerical analysis we take gy = gy> (the 
U(l)' coupling constant), which does not conflict with the unification of the gauge couplings. 

From the tree-level potential one can derive the minimization equations for the VEVs 
v u , Vd, Vs and the phase 9s(9 s ). These relations yield conditions relating the VEVs to the 
physical Higgs masses. 



The spectrum of physical Higgs bosons consist of three neutral scalars (h, H, H'), one CP 
odd pseudoscalar (A ) and a pair of charged Higgs bosons if ± in the CP conserving case. In 
total, the spectrum differs from that of the MSSM by one extra CP-even scalar. Notice that, 
the composition, the mass and the couplings of the lightest Higgs boson of U(l)' models 
can exhibit significant differences from the MSSM, and this could be an important source of 
distinguishing signatures in the forthcoming experiments. It is important to emphasize that 
these models can predict naturally larger values for nth, the lightest neutral Higgs boson 
masses, which are more likely to agree with the boson mass seen at the LHC. While we can 
safely require > 90 GeV for all numerical estimates [24], in principle, it is possible to 
obtain larger values such as ~ 140 GeV within some of the E 6 based models. In our 
evaluations, we shall impose rrih ~ 124 — 126 GeV, in agreement with the mass of the particle 
observed at the LHC. 



B. One-loop Corrections to the Higgs Potential 



The tree level potential in Eq. is insufficient to make precise predictions for masses 
and mixings, and thus we include loop corrections. For this we use the effective potential 
approach. Not all of the CP violating parameters are free parameters, and loop corrections 
induce certain relationships among them. The one-loop corrected potential has the form 
V = Vtrp.p + AV, where y ree is defined in ([6l), and Ay is the one-loop Coleman- Weinberg 



potential 25] : 

Ay 



64vr 2 



Zj(-l) 2J+1 (2J+l)M*(H u ,H d ,S) 



In 



M 2 (H u ,H d ,S) 3 
A 2 2 



(8) 



where hA represent the mass matrices of all the particles in the theory. While many particles 
and their superpartners could be added for the calculation of the loop corrections, we include 
here the dominant contributions coming from the top and bottom sectors (/ = t, b) for both 
the quarks and scalar quarks, so that both contributions from small and large tan /3 values 
can be investigated safely. Specifically 



Ay 



6 



64vr 2 



EE 



• 2 \2 



In 



2 

m~ 

Jk 

A 2 



2(m 



2\2 
f) 



( rrtf 



(9) 



f=b,t Kk=l,2 

In this expression the masses depend explicitly on the Higgs field components: for instance 
the bottom mass-squared is given by m 2 = Y?\H d \ 2 , and the top by m 



Y t 2 \H u \ 2 , and 
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the scalar quark masses-squared are obtained by diagonalizing the mass-squared matrix, the 
unitary matrix Sf as SjM 2 Sf = diag(m~ , m~), with f — t,b. 

The vacuum state is obtained by requiring the vanishing of all tadpoles and positivity of 
the resulting Higgs boson masses. The vanishing of tadpoles for V along the CP-even direc- 
tions (pH u ,H d ,s an d CP-odd directions (p Utd ,s allows the soft masses m 2 HuHdS to be expressed 
in terms of the other parameters of the potential. The tadpole terms are obtained from 

*-(».■ 

where means that we evaluate the derivative at the minimum of the potential, V = 
V tree + AV , and $j = <p u , <p d , <p s , <p u , <£d, <Ps- Since all tadpole terms must vanish, enforcement 
of 7i,2,3 — is used to obtain rriH u ,H d ,s, respectively, and 74,5,6 can be used for the phase of 
the trilinear coupling (As), which is 6 S - In fact at the tree-level the result is 9 S = 0, but 
loop corrections induce this quantity to be non-zero. For instance, at the tree level, using 

7i, 7*2 and % (given explicitly in the Appendices), one can express Higgs mass-squared as 
, AsY s cos(6z + 6 s )v d v s Q Hu U + Y 2 (v 2 + v 2 ) g^g-tj) 

2 AsY s cos(9s + 9 s )v u v s Q Hd Yi + Y[ (vj + vj) g 2 (v 2 u -vj) 
mHd ~ V^ d 2 8 ' (12) 

2 _ A s Y s cos(6x + 6 s )v d v u Q S U + Y 2 (v 2 + v 2 u ) 

ms ~ T^s 2 ' { } 

where 

n = qUQhA + Q s v 2 s + Q H y u ). (14) 

At tree-level %, % and % are zero, but at one-loop level they all induce the same non-zero 
result. We collected the full form of the tadpoles 71, % and % in the Appendices. Using the 
tadpoles along the CP odd directions, the phase of the trilinear coupling of S (As) emerges 
as a radiatively induced quantity, 

9 -> ,, n -i( ^ F ^A b Y 2 + F t S t A t Y 2 ) \ 

9s ^-sm ^ J-0 S , (15) 

where we defined S t = sm(9 t + #e) and S& = sin(# 6 + 9%). We define cosine of the same 
quantities: C t = cos(9 t + 0%) and C b = cos (9 b + 9%). Here F t and F b are loop functions: 



2 2 

m r m 



F / = -2 + .„|^j-l„(-|jfi, (16) 

where / = t,b refers to top and bottoms and we defined Q as the SUSY breaking scale, 
A f — m 2 r - m 2 r and E / = m| + m 2 ~ t . 

J J2 h J 72 Jl 
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C. The Higgs Mass Calculation 



We now turn to the Higgs mass calculation at one-loop in the presence of CP violation 
in the stop and sbottom LR mixing. The mass-squared matrix of the Higgs scalars is 

M 2 - = ( ^ v) . (17) 

In the above $; = (0j,<^j). Two linearly independent combinations of the pseudoscalar 
components tp U) d,s are the Goldstone bosons Gz and Gz>, which are used to give mass 
to the Z and Z' gauge bosons, leaving one physical pseudoscalar Higgs state A , which 
mixes with the neutral Higgs mass states in the presence of CP violation. In the basis of 
scalars B = {0 U , <fi d , <p s , A }, the neutral Higgs mass-squared matrix M 2 takes the following 
symmetric form 

/ M 2 U M\ 2 M 2 3 M\, \ 



M 2 



HO 



M\ 2 M 2 2 M 2 23 M 2 , 
M 2 3 M 2 23 M 2 3 M 2 34 

<2 A^2 A^2 AW2 



(18) 



\ M{, M 2 U Mi, M\, ) 
The mass-squared matrix can be diagonalized by a 4 x 4 orthonormal matrix O. In doing 
this we follow the convention OM 2 H oO^=di&g(m 2 H o,m 2 H o,m 2 H o,m 2 H o), where, to avoid dis- 
continuities in the eigenvalues, we adopt the ordering: m H o < m H o < m H o < m H o. The 
elements of O determine the couplings of Higgs bosons to the MSSM fermions, scalars, and 
gauge bosons. 

The results for the entries of the neutral Higgs (mass) 2 matrix are collected in the Ap- 
pendices. As an example, we show here one of the masses for the CP-conserving case. When 
CP is conserved all M 2 , and M\ entries should vanish, with the exception of the M 2 ,, term, 
which is actually the pseudoscalar Higgs (mass) 2 term. When CP is conserved M\ is 

VdVgV u v d v%v u 
where co 2 = v 2 v 2 s + v 2 d v 2 u and re = 3/(32tt 2 A 2 A^). 

Calculation of masses of the charged Higgs bosons is very similar to the neutral ones and 
we obtain the following mass-squared matrix 

/ M\t M\t » 

M 2 H± =\ 11 12 , 20 

1 M\t M 2 2 ? 1 
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and the eigenvalue of this matrix yields, when CP is not conserved, the expression 



w\ 2 A 2 

= o 2^ 6 2V (M3Y b 2 vUF b X b (»A b (C b (v 4 d + v 4 u ) + 2S b v 2 v 2 u ) - A\v d v\ - ^v 3 d v u ) 
- Y> 2 b v d vl{F b + G h - 2) + A 2 b (G b - 2)v d v 3 u ) - ^(v* + i£)(87rV>«(V - g 2 2 v 2 s ) - fi X 4) 

9 

in 

+ GY b ^ b vy s vl(\n(^) - 1)) + 3Z b Y?vl(F t Z t (»A t (C t (vj + v'J + 2v 2 S t vl) - A\v\v u - fv d v\) 



Q 

Q 



2 

mi 



- vl\Y? t v u (F t + G t -2) + v 3 d (G t - 2)A 2 t v u ) + 6S^ d 3 F>lE^(i n (_A) _ 1)). ( 2 l) 



where we defined the loop function 




G / = 2 + ln -f \f±, (22) 



with f — t,b. From this it is easy to obtain the mass of the charged Higgs in the CP 
conserving case. This can be achieved by taking the limits C t — > 1, C b — > 1 and St — > 0, 
S b — > 0. We present explicitly the four entries of the charged Higgs mass-squared matrix in 
the Appendices. 



D. The Neutralino Mass Matrix in U(l)' 

The presence of the CP-violating affects the chargino, neutralino and scalar quark mass 
matrices. As we are concerned here with the (tree-level) Higgs decays into neutralinos, we 
show the effect on the phases on the neutralino mass matrix. Note that the chargino mass 
matrix is unchanged from the MSSM one, though it depends on U(l) ; breaking scale through 
the fi — > fi e ff parameter in the mass matrix. Similarly, the elements in the sfermion mass 
matrices are modified due to the presence of the Z' boson. Their explicit expressions have 
appeared elsewhere 26]. 



The neutralino sector of the U(l)' is like the MSSM, but enlarged by a pair of higgsino 
and gaugino states, namely S (referred to as singlino) and B', the bare state of which we 
call bino-prime, while Z' (zino-prime) is the physical mixed state. The mass matrix for the 

11 



six neutralinos in the (B, W 3 , H%, H®, S, B') basis is given by a complex symmetric matrix: 



M l -M z c p s w M z s p s w M K \ 

M 2 M z cpc w -M z s p c w 

-M z c p s w M z cpc w -Mf -VxSp Q,H d M v c p 

M z s p s w -M z s p c w -faff ~H\Cp Q Hu M v Sp 

-fixsp -Hxcp Q S M S 

M K Q Hi M v c p Q Hu M v s p Q S M S M[ ) 



(23) 



with gaugino mass parameters M 1 , M 2 , M[ and [27] for B , W 73 , B' and B — B' 
mixing respectively, tan/3 = v u /vd, and 0^ denotes the electroweak mixing angle. After 
electroweak breaking there are two additional mixing parameters: 



M v = qy'V and M s = gY'Vs- 



(24) 



Moreover, the doublet-doublet higgsino and doublet-singlet higgsino mixing mass mixings 
are generated to be 



a ff — Y*—e i6s 
y2 



V2 ' 



(25) 



where v = yv%~+ v\. The neutralinos mass eigenstates are Majorana spinors, and they can 
be obtained by diagonalization 



X° = Mjipj , Xo = (Xo, X°i) T 



(26) 



The neutralino mass matrix is diagonalized by the same unitary matrix 



J\f'M x oJ\f = diag(m x o, rh x o] 



(27) 



The additional neutralino mass eigenstates due to new higgsino and gaugino fields encode 
the effects of U(l)' models, wherever neutralinos play a role such as in magnetic and electric 
dipole moments, kaon mixing, or in Higgs decays. 
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III. CONSTRAINTS AND IMPLICATIONS FOR THE CP VIOLATING HIGGS 
SECTOR 



A. Electric Dipole Moments 



The experimental bounds on the electric dipo 



cm and the electron d e < 1.8 x 10 27 e cm 17|, |29j, are some of the most tightly bound 



e moments of the neutron gL < 6.3 x 10 26 e 



measurements in physics. The electric dipole moment (EDM) of a spin-| particle is defined 



from the effective Lagrangian 



28 



d = --drfa^F'", (28) 



and it is induced at the loop level if the theory contains a source of CP violation at the 
tree level. Unlike the SM, where the EDMs are generated through the phase of the CKM 
matrix at higher loop level and are thus small, in MSSM, where they are generated at 
one- loop level, the electric dipole moments are very important, and they provide important 
restrictions on the parameter space of the model. In U(l)' supersymmetric models, they 
acquire contributions from gluinos (for neutron EDM) and chargino and neutralino (for 
both neutron and electron EDMs), and the contributions are generated by /i e // = fie t0s , 
with an additional contribution generated by the Z' neutralino. The EDM was analyzed in 



22l | in the limit in which the sfermions are much heavier than the Z '. 

The neutralino contributions to EDMs tends to be overall subdominant. To suppress 
the EDMs we can proceed as in the MSSM [28]: we can require that the trilinear stop 
coupling be mostly diagonal A l t ~^ 3> A 1 ^ (that would suppress the sfermion contribution); 
we can assume cancellation between different SUSY contributions (in particular destructive 
interference between gluinos and charginos); or we can require that the first and second 
generation sfermion masses be in the TeV region. Alternatively, one can assume generically 
small CP-violating phases, a path we do not wish to follow here, not just based on natural- 
ness, but because we wish to investigate the effects of the phases on Higgs phenomenology. 
In the case where gyi = gy, the case we consider here, the constraints on U(l)' parameters 
are similar to those on the MSSM. The parameter space we choose for our benchmark points 
insures that the contributions to the EDMs are sufficiently small. 
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B. CP violation in K° - K° mixing 

The physical phases of the Higgs singlet and in the scalar fermion, chargino and neutralino 
mass matrices could alter the the value of the measure of the CP violation in K° — K° mixing, 



measured to be e K = (2.271 ± 0.017) x 10~ 3 (l7|. 



The contributions to the indirect CP violation parameter of the kaon sector, defined as 

with Am, the long- and short-lived kaon mass difference, and M.yi the off-diagonal element 
of the neutral kaon mass matrix, is related to the effective Hamiltonian that governs AS" = 2 
transitions as 

/ TfO\njAS=2\ ff0\ 

M * = 2m ' With ^e A / =2 = E c ^- ( 30 ) 

i 

Here q are the Wilson coefficients and Oi the corresponding four-fermion operators. In the 
presence of SUSY contributions, the Wilson coefficients can be decomposed as a sum 

where the first contribution is the SM one, the second is the charged Higgs, and the rest are 
supersymmetric contributions. In U(l)' models, the dominant supersymmetric contributions 
come from the chargino mediated box diagrams, and the AS" = 2 transition is largely 
dominated by the (V — A) operator 0\ = dj^Pisdj^PLS, similar to the MSSM, and the 
chargino contribution is larger than the charged Higgs contributions. The contribution in 
terms of the bare chargino states is approximately [23|: 

x |A^sin^ (m|) la I{ r w±^ r H±^t L ,n R )} , (31) 

where is the kaon decay constant and tjik the kaon mass; Vij are the Vckm elements, 
(rriq) is the average squark mass, taken equal to Msusy; m w± — M 2 is the wino mass, and 
mg± = /i is the higgsino mass, and rj = mf/(mq) 2 . The non-universality in the LL soft 
breaking masses is parametrized by (^|) 12 , and the non-universality in the soft trilinear 
terms is parametrized by AA t = A] 3 — A 23 . Finally, I is the loop function which can be 
reduced to elementary functions in the limit of degenerate squark masses. 
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Scanning the parameter space of the model, we checked that one can find parameter sets 
that satisfy the minimization conditions of the Higgs potential, have an associated Higgs 
boson spectrum compatible with the LHC boson, and still succeed in obeying the bound 
for the observed value of Ek- From Eq. (1311) . it appears that Ek depends on 1/M| USY . To 
satisfy the experimental value of Ex, values of Msusy > 1 TeV would have to be assumed, 
or A A t = A] 3 — Af 3 «C 1, in agreement with our EDM considerations. Too small values of 
Msusy might generate a light Higgs boson spectrum already excluded by LEP and LHC, 
and for Msusy ~ 1 TeV, which is consistent with our squark and slepton masses, the 
supersymmetric contributions to Ek are consistent with the experimental constraints. 



IV. NUMERICAL ANALYSIS 



As mentioned in the introduction, gauge extensions of the SM by one or several non- 
anomalous U(l)' gauge groups can arise naturally from a string-inspired EqSSM model 
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lo| . In EqSSM models the matter sector includes a 27-representation for each family 
of quarks and leptons (including right-handed neutrinos), Higgs representations (doublets 
H u and and singlet S), and three-families of extra down-like color triplets. Anomaly 
cancellation occurs generation by generation, and gauge coupling unification requires another 
pair of Higgs-like multiplets. Breaking of E e yields SU(3) x SU(2) x U(l)y x U(l)' as a 
low energy group. Anomaly-free 11(1)' groups are thus generated this way, directly, or as 
a specific linear combination. We first define the models that shall be investigated in our 
numerical analysis. They al 
anomaly free groups U(l)^ 



emerge from breaking of higher groups [31|. For instance, the 



32] and U(l) x [33] are defined by: 



E 6 ->• SO (10) x U(l)^, SO (10) ->• SU(5) x U(l) x . 

In general a U(l)' = U(1)e 6 group is defined as U(1)e 6 = cos9e 6 U(1) x + sm9 Efi U(l)^, and 
we distinguish among the different scenario by the values of 9e 6 '- 

9 V = n — arctan w | for U(l) v wich occurs in Calabi-Yau compactification of heterotic 



strings 



34|; 



9s = arctan y/TE/9 for the secluded U(l)s, where the tension between the electroweak 
scale and developing a large enough Z' mass is resolved by the inclusion of additional 



singlets 



ad; 
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9 j = arctan J | for the inert U(l)j, which has a charge orthogonal to Q v [36]; 



masses 



7T 
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On = arctan y 15 for U(1)n, where v c has zero charge, allowing for large Majorana 



and 



• Oif> — ~n for U(l)^, defined above from the breaking of E§ |32|. 

In the Table [J below, we list the charges for the fundamental representations of Eg in the 
U(l)' models which we use for numerical investigations of Higgs boson properties. 



SO(10) representations 


SU(5) representations 


2v / 15Qr ) 


2Vl5Qs 


2Qj 


2VWQ N 


2VQQ^ 


16 


10 (u,d,u c ,e + ) 


-2 


-1/2 





1 


1 


(u, d, v, e~ , u c , d c , v c , e + ) 


5* (d c ,u,e-) 


1 


4 


-1 


2 


1 




v c 


-5 


-5 


1 





1 


10 


5 {H u ) 


4 


1 





-2 


-2 


(H u , H d ) 


5* (H d ) 


1 


-7/2 


1 


-3 


-2 


1(S) 


1 (S) 


-5 


5/2 


-1 


5 


4 



TABLE I. Values of f7(l)^, U(l) s , U(X)l, U(1) N and U(l)^ charges for the 27 fundamental rep- 
resentation of Eq decomposition under SO(10) and SU(5) representations. The charge for each 
model is defined as Q = cos 9 e 6 Qx + sin Oe 6 QiI>- 

In what follows, we investigate the consequences of each of the anomaly-free groups 
on the Higgs production and decay at the LHC In Table HT1 we list the relevant benchmark 
parameters for each of the choices, for both the CP-violating (CP-conserving) Higgs sector J3- 
In addition to the phase 8 S (which defines the CP violating scenario of each model), the 
values of tan/3 and of fi, we give the U(l)y and SU(2)l gaugino masses Mi and M2, the 
left and right handed squark soft mass parameters Mq 1 and Mjj i (all taken to be 1 TeV, 
including the masses in the down scalar sector, not explicitly shown), the trilinear couplings 
in the top and bottom scalar quark sectors, A t and A&, and the ratios Ry> = -r~r an d 
M 

Ryy = ~tj~i as defined in [26]. The constraints on the mass parameters, constraining the 
choice of benchmark values, are: 



By CP violating scenario, we mean the specific case where 9 S is given by the values in Table [TTJ 
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Requiring the lightest Higgs mass to be very close to 126 GeV, in agreement with the 
ATLAS and CMS results; 

Requiring the next lightest neutral Higgs boson to have mass m H o > 600 GeV (as it 
has not been observed at LHC); 

Requiring the lightest neutralino mass to be consistent with collider limits on Z boson 
decays, but also to allow for the possibility of the neutral Higgs boson to decay into a 
neutralino pair; 

Choosing the lightest neutralino to be the LSP and requiring that the relic density 
constraint be satisfied; 



Choosing the Z' boson mass to be consistent with present limits 17j; 

Choosing scalar masses and trilinear couplings which satisfy constraints from EDMs 
and CP violation in the kaon sector, as described in the previous section (IHIj) . 



As we would like to allow the Higgs boson to be kinematically allowed to decay into two 

TiZ -> xlxl) < 3 



43 



neutralinos, we impose the LEP constraint on the Z boson width 
MeV. This constraint allows for a weakening of the Particle Data bound |17| . especially as we 
do not impose the supersymmetric grand unified theory relationship M\ = (5/3) tan 2 Qy^Mi, 
and allow M\ and M2 to be free parameters, as given in Table [TXJ Note that in particular, 
the bino mass is chosen to be light to allow Higgs decays into neutralinos, while value for 
M 2 insures that the chargino mass will be m~± > m H o/2. The scalar fermions are heavy 
to satisfy bounds from the EDMs and 6k- We choose the value of 9 S for each model to 
maximize the invisible decay width for the lightest Higgs boson, while satisfying the other 
constraint^. 

Based on the input parameters, we calculate the spectrum of the physical masses of the 
extra particles in the model, which are used in our numerical evaluations. These values are 
given in Table IIHI We also included in this table the relic density of the dark matter for 
all scenarios. Throughout our considerations the lightest neutralino Xi i s t ne lightest super- 
symmetric particle (LSP) and thus subject to cosmological constraints. The relic calculation 



* Our benchmarks are different from those NMSSM 1JJ. where CP conservation was assumed, and where 
the dominant decay mode of the lightest CP-even Higgs is into the pseudoscalar Higgs boson pairs. 
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Parameters 


u(i) v 


U(l)s 


U(l)i 


U(1) N 




O s 


42(0) 


75(0) 


60(0) 


55(0) 


33(0) 


tan (3 


1.8(1.7) 


1.46(1.42) 


1.3(2.5) 


1.5(1.8) 


1.75(1.75) 




360(360) 


715(730) 


465(461) 


292(295) 


285(290) 


Mi 


48(50) 


56(59) 


57(50) 


49(50) 


49(51) 


M 2 


125(130) 


115(120) 


135(120) 


130(170) 


140(160) 


MQi 


1000(1000) 


1250(850) 


750(600) 


2000(300) 


1000(1000) 


M Q2 


1000(1000) 


1250(850) 


750(600) 


2000(300) 


1000(1000) 


M Qz 


1000(1000) 


1250(850) 


750(600) 


2000(300) 


1000(1000) 


M Vl 


1000(1000) 


1250(850) 


750(600) 


2000(300) 


1000(1000) 


M U2 


1000(1000) 


1250(850) 


750(600) 


2000(300) 


1000(1000) 


M U3 


1000(1000) 


1250(850) 


750(600) 


2000(300) 


1000(1000) 


\M 


1850(2000) 


2200(2500) 


2500(1500) 


2250(2000) 


2000(2000) 


\A b \ 


2000(2000) 


2500(2500) 


2500(1500) 


2500(2000) 


2000(2000) 


R Y > 


1(1) 


2.5(0.1) 


0.1(6.6) 


1(1) 


5(5) 


Ryv 


1(2.2) 


0.1(0.1) 


2(6.6) 


6(2.7) 


0.1(5) 



TABLE II. The benchmark points (in GeV) for the CP-violating (CP-conserving) 11(1)^, U(l)s, 
U(l)j, U(1)n and t/(l)i/> versions ofU(l)' models. 



is straightforward using the Micromegas package [39], once we include the U(l)' model files 



fromCalcHEP 40j. All the numbers are within the lcr range of the WMAP result 
the Sloan Digital Sky Survey 



42 1 



4lj from 



0.111 



-0.011 
-0.015 



(32) 



The relic density of the dark matter f2rjM^ 2 i s very sensitive to the free parameter Ry> listed 
in Table [Til As the lightest neutralino plays an essential role in the decay of the lightest 
Higgs boson, we first show the dependence of its mass, and of the relic density with the 
CP violating parameter 9 S in Fig. [TJ In all of the U(l)' models under study the lightest 
neutralino is mostly bino. The variations of its mass and of the relic density with the other 
CP violating phase 9 t are negligible. Note that the mass of the LSP increases smoothly 
with increasing 8 S , while the relic density measurement (shown as a green band in the right- 
handed part of the plot) poses restrictions on the combined LSP mass and CP violating 
parameter. The values of 9 S for various models listed in Table [TTJ fall into the range of the 
values allowed by the relic density (within the green band). We incorporate these restrictions 
in our analysis of Higgs mass and decay widths. 
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Masses 


U( 


l)r? 


U(l)s 


U( 


1)/ 


U(1) N 


U(l)i> 


m Z ' 


1510 


4510) 


1507(1539) 


1513 


1500) 


1502 


1517) 


1513(1540) 


Xi 


43 


;43) 


55(55) 


54 


14) 


43 


12) 


42(42) 




108 


109) 


112(110) 


125 


107) 


111 


137) 


114(128) 


A3 


361 


[361) 


715(730) 


464 


163) 


292 


'297) 


286(292) 




386 


;388) 


726(742) 


485 


179) 


326 


'336) 


322(331) 


A5 


1487 


1489) 


1440(1536) 


1514 


1378) 


1505 


1498) 


1396(1438) 




1535 


1540) 


1580(1543) 


1521 


1711) 


1 556 

laUU 




1641(1694) 


m~± 

Xi 


107 


107) 


111(110) 


124 


106) 


108 


134) 


111(125) 


m-± 


382 


;384) 


724(740) 


481 


177) 


321 


'332) 


318(326) 


Vfl ttO 
n l 


125.0 


125.0) 


125.6(125.0) 


125.8 


126.0) 


125.6 


126.0) 


125.4(125.0) 


n 2 


743 


747) 


969(1027) 


788 


;930) 


642 


'688) 


665(679) 


TUttO 

n 3 


750 


754) 


977(1033) 


798 


;933) 


652 


'695) 


673(687) 


"4 


1510 


1510) 


1508(1539) 


1513 


1500) 


1502 


1517) 


1513(1540) 


m H ± 


572 


;543) 


717(711) 


507 


;802) 


418 


'504) 


486(486) 


mg L 


1341 


1341) 


1837(1616) 


1306 


1219) 


700 


742) 


1134(1139) 


rrie R 


1054 


1054) 


1154(695) 


748 


;598) 


513 


'564) 


1133(1137) 


rrinr 


1341 


1341) 


1837(1616) 


1306 


1219) 


700 


742) 


1134(1139) 


m p.R 


1054 


1054) 


1154(695) 


748 


;598) 


51 3 


'564) 


1133(1137) 




1054 


1054) 


1154(695) 


748 


;598) 


513 


'564) 


1133(1137) 


m f2 


1342 


1341) 


1837(1616) 


1306 


1219) 


700 


742) 


1135(1139) 




1340 


1340) 


1836(1615) 


1306 


1217) 


699 


739) 


1133(1137) 




1340 


1340) 


1836(1615) 


1306 


1217) 


6QQ 


73Q) 


1133(1137) 




1340 


1340) 


1836(1615) 


1306 


1217) 


699 


739) 


1133(1137) 




1054 


1054) 


879(874) 


999 


;998) 


1106 


1108) 


1133(1137) 


m u R 


1055 


1055) 


882(877) 


1001 


1000) 


1107 


1109) 


1134(1138) 


&L 


1056 


1055) 


880(875) 


1000 


1001) 


1107 


1109) 


1134(1139) 


a R 


1340 


1340) 


1675(1698) 


1463 


1457) 


1203 


1207) 


1133(1138) 




1054 


1054) 


879(874) 


999 


;998) 


1106 


1108) 


1133(1137) 


m- CR 


1055 


1055) 


882(877) 


1001 


1000) 


1107 


1109) 


1134(1138) 




1056 


1055) 


880(875) 


1000 


1001) 


1107 


1109) 


1134(1139) 


m SR 


1340 


1340) 


1675(1698) 


1463 


1457) 


1203 


1207) 


1133(1138) 




919 


;9ii) 


659(670) 


788 


;894) 


938 


;968) 


994(1002) 




1201 


1207) 


1085(1070) 


1200 


1122) 


1277 


1275) 


1281(1283) 


mr 

01 


1056 


1055) 


880(875) 


1000 


1001) 


1107 


1109) 


1130(1135) 


mr 


1340 


1340) 


1675(1698) 


1463 


1457) 


1203 


1207) 


1137(1141) 




0.114 


;0.120) 


0.100(0.102) 


0.113 


;o.i2o) 


0.111 


;o.ii7) 


0.117(0.101) 



TABLE III. The mass spectra (in GeV) and the relic density Odm values for the CP-violating 
(CP-conserving) version of the scenarios considered given in Table [7I| for the U(l)' models. 
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FIG. 1. Mass of the lightest neutralino and the relic density as functions of 9 S (the phase of the 
new singlet S) for the CP-violating versions of "J7(l)^, U(l)s, U(V)i, 17(1) jv and U(l)^ models. The 
green band indicates the experimentally allowed region. 

We proceed to examine the effects of the CP violating phases on the masses, production 
cross sections and branching ratios of the lightest Higgs boson. 

A. The lightest CP-even neutral Higgs boson 

The observation of the new boson at the LHC has fueled speculations of its nature (is 
it or not the SM Higgs boson), coupled with analyses of its mass and couplings, and their 
comparison with the experimental data. ATLAS [45] and CMS [46] have reported updates 
on the combined strength values for main channels, including H° — > bb, 77, r + r~, WW*(—> 
Ivlv) and ZZ*(—> 4£). While the results still have significant experimental and systematic 
uncertainties, these are expected to decrease with LHC operating at y/s = 14 TeV and 
increased luminosity. The precise determination of the Higgs couplings to different channels 
will establish whether the boson observed at the LHC is the SM Higgs boson. In our 
analysis, we wish to explore the possibility that Higgs boson decays in a non-SM fashion, in 
particular, that it can decay significantly invisibly. A invisible decay mode is very hard to 
measure directly at the colliders. However, it is not difficult to be inferred indirectly. The 
total decay width of a SM Higgs boson with mass of 125 GeV is approximately T^o = 4.2 
GeV. A discrepancy between the theoretical and experimental value for the width would be 
an indication of additional decay channels beyond SM. Similarly, reduced decay branching 

20 



80 100 
S 



u N 

120 140 160 180 



ratios into known SM Higgs decay modes, in particular for H° — » bb and H° — > t + t~ 
(dominant for m#o = 126 GeV), could also indicate that other decays are important. 

At tan (5 ~ 1 the lightest Higgs mass is determined mostly by the new F- and D-terms 
in the Higgs tree-level potential, and is thus sensitive to the trilinear Yukawa coupling 
and the gauge coupling gY>{= 9y) i n the numerical analysis. We first present our results for 
the dependence of the masses on the CP violating phases arg (/i e //) = 9 S and arg (A t ) = 9 t , 
as well as with tan j3 in Fig. [2j One can see that the mass variations with 9 S and 9t are 
significant, especially in U(l)s, where large regions of the parameter space for both phases, 
if combined with other measurements, can be eliminated. The dependence on tan (3 from 
the third panel of the figure seems to indicate that only low values tan w 1 — 2 are allowed 
for all U(l)' models, in agreement with the values chosen in Table [Hj 
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FIG. 2. Mass of the lightest neutral Higgs boson as a function of 9 S (the phase of the new singlet 
S), 6t (the phase of the soft coupling At), and tan/3 for the CP-violating versions of 11(1)^, U(l)s, 
U(1) N and 17(1)^ models. 

To analyze the decay width of the lightest Higgs boson, we first calculate total production 
cross section of the lightest Higgs boson (H®) in various models in Table IIVt for 9 S = (no 
CP violation) and for 9 S as in Table HT1 (with CP violation). We list associated Higgs- vector 
boson cross sections, and the total cross section for the vector boson fusion. Though sub- 
dominant production modes for Higgs bosons, these are the dominant channels for observing 
an invisible decay of the Higgs boson [ 1 1 j] . Note that we do not include here the dominant 
production mechanism gg — > H®, as this mode is plagued by large QCD corrections, and 
thus it is difficult to isolate the invisible decay of the Higgs boson, which in this production 
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channel is expected to come from gg — > H^+iet, and be small. As expected, the vector-boson 
fusion production mechanism dominates over the Higgs-vector boson associated production 
in all models. The numbers are fairly consistent across the models, and largely independent 
of CP violating phases. Thus we forgo plots of the production cross sections and expect that 
any differences would show up in the branching ratios of the lightest Higgs boson. 



Observables 




U(l)s 


U(l)i 


U(l) N 




a(pp -> H^Z) 


639(642) 


631(647) 


628(610) 


628(624) 


634(642) 


a(pp -> H?W + ) 


720(725) 


708(725) 


705(687) 


708(701) 


711(720) 


a(pp^H?W-) 


445(447) 


437(448) 


435(424) 


437(433) 


439(444) 


a(pp^H?jj(VBF)) 


4983(4930) 


4848(4920) 


4861(4840) 


4874(4850) 


4873(4893) 



TABLE IV. Total cross sections of associated production channel (H®X ) and vector boson fusion 
production channel (H®jj) (in fb) for the CP-violating (CP-conserving) versions of [7(1)„, U(l)s, 
U(l)i, U(1)n an d U(1)j, models considered in the paper. 



We list the dominant decay branching ratios (in %) for the lightest neutral Higgs in 
our model and for comparison, in the SM in Table |V] again for no CP violation (9 S = 0) 
and with CP violation (with phases as given in Table [TTJ) . The branching ratios, as well 
as the cross sections are largely independent of the 9 t phase. One can see that, while the 
production cross sections are fairly independent of the CP violating phase 9 S , the branching 
ratios are not, showing significant differences between the various U(l)' scenarios and the 
SM in the branching ratios. First, given the fact that the lightest neutralino (the LSP) has 
mass m^o < m H o/2, the Higgs boson has a considerable branching ratio into XiXii that is, 
a significant invisible width_j|. This is accompanied by a reduction in the branching ratio to 
other two-body decays, in particular r + r~ and WW*. Of all the U(l)' scenarios, the invisible 
width is the smallest in U(l)s, though comparable with the decay width into t + t~ for the 
case of no CP violation. For the other U(\)' models, the branching ratio for the invisible 
decay goes from a low 9% in U(l)i with no CP violating phases, to 54% in U(l)^ with CP 
violation. A general feature emerging from Table |V] is that the invisible width is enhanced 
in the presence of CP violation [9 S ^ 0) over the case with 9 S = 0. This is particularly 
§ Note that in principle the Higgs boson can decay into sneutrinos, which can then cascade into neutralinos, 
contributing to the invisible width. We preclude this possibility here, as m„ < m H o/2 would require soft 
left- handed slepton masses of 0(100) GeV, in conflict with the EDM constraints. 
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strong in the case of U(l)s, where the branching ratio increases for 9 S (as in Table ITT]) to 3 
times of its CP-conserving value, and for U(l)j where it increases more than twofold. The 



decay into t 
the MSSM 



re invisible mode can reach over 50%, which is similar to the value obtained in 
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Note that the decay into invisible modes is sometimes at the expense of 
the main SM decay into bb. In two of the models studied, U(l)s and U(l)j the H® — > bb 
branching ratio is in fact increased with respect to the SM value, while in U(l) v , U(1)n and 
U(l)i/j it is suppressed with respect the SM expectations. But a general feature of all these 
models is the strong suppression of the H® — > (W + W*~ + W* + W~) and H — > t + t* decay 
modes, expected to have a branching ratio of 21.5% and 6%, respectively, in the SM, but 
much smaller here. The branching ratio for the decay H® — > t + t~ is between ~ 2 — 3.5%, 
while that for — > WW* ranges between ~ 5.5 — 12%. In a nutshell, the Higgs decay 
into the invisible mode X1X11 is & t the expense of — > W + W~ and t + t~ in all models, 
and occasionally due to a suppression of — > bb in some models. This behavior is not 
unexpected, as previous studies have indicated that for light Higgs masses, the decay into 
neutralinos and Higgs pseudoscalar pairs (if kinematically allowed) dominate, at the expense 
of the SM decay modes. Increasing the lightest Higgs mass opens allowed channels, but the 
branching ratios are affected by the mixing with the singlet Higgs field, the pseudoscalar 
and the effect of the CP violating phase. However, due to differences in decay patterns 
among various anomaly-free versions of the U(l)' models, a more precise measurement of 
the Higgs boson branching ratios at the LHC will serve not only to differentiate between 
the SM and the U(l)' model, but among the different versions of U(l) n s. In Fig. [3] we 



Branching 


I Ratio 




U(l)s 


U(l)i 


U{1) N 




SM 


BR(Hj ) -S 


• xm) 


36.0(34.0) 


8.0(2.6) 


20.0(9.0) 


49.0(41.0) 


54.0(42) 




BR(H5 - 


bb) 


48.0(49.0) 


70.0(73.0) 


60.0(66.0) 


38.0(44.0) 


36.0(43.0) 


60 


BR(H? 


T~T+) 


2.3(2.4) 


3.5(3.6) 


3.0(3.3) 


1.9(2.2) 


1.8(2.2) 


6 


BR(H? 


WW*) 


7.4(7.2) 


10.9(11.1) 


9.8(12.0) 


6.1(7.5) 


5.3(6.6) 


21.5 



TABLE V. Dominant branching ratios (in %) of decay channels for the CP-violating (CP- 
conserving) version of the 1/(1)^, U(l)s, U(l)i, U(1)n and U(l)^ scenarios considered, and in the 
SM. 



plot the variation of the branching ratios of the lightest Higgs boson with the CP violating 
phase 6 S . In the first two panels, we depict the dependence of the BR(f/^ — > XiXi) with 9 S 
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and tan/3. As we have seen previously tan/3 ~ 1 — 2 (as in Table |TT]), and in that region 
the invisible decay width is large, and very sensitive to tan 0. We show the variation of 
the branching ratios of the other dominant SM and U(l)' modes, as well as the that for 
BR(H® — > A£), because the LHC is sensitive to this decay in the 124-126 GeV mass range, 
and the value is expected to become more precise. Note that we did not include any of 
the loop dominated decays, such as H® — > gg 1 r j r j ) as these are sensitive to the masses and 
mixing parameters of the (numerous) particles in the loop, and there no new contributions 
to these processes with respect to MSSM. The third panel in the top row of the figure shows 
that, while the BR(H® — > bb) in SM seems to fall somewhere in the middle of predictions 
for U(l)' models, the SM BR(/f° — > r + r~) (bottom row, left side) is 6%, and outside the 
range of U(l)' models, and so is the SM value for the BR(H® — > WW*) (bottom row, middle 
panel). BR(f/^ — > 4£) (bottom row, right panel) is also beyond the upper high end of U(l)' 
models predictions; the value expected in the SM is 0.013% while in the U(l)' models, the 
BR's fall in the ~ 0.0025 — 0.0055% range. The results for these decay widths might be more 
meaningful experimentally than the invisible Higgs width, which is difficult to measure. 

B. The second lightest neutral Higgs boson 

If the underlying symmetry in nature is not the SM, it is very likely that more Higgs 
boson states will be observed. The U(l)' models all predict additional neutral and charged 
Higgs states. The present collider bounds indicate that the mass of the second lightest 
Higgs boson must be heavier than about 600 GeV. In our model, this mass shows explicit 
dependence on the CP violating phases 8 S and 8 t . This dependence is correlated with the 
lightest boson mass. As the m# > 600 GeV mass region will be available to LHC working at 
increased -y/i = 14 TeV, we show the mass dependence of the second lightest neutral Higgs 
boson in Fig. HI The variation of this mass with either on the CP violating phases 9 S or Q t 
is not as pronounced as for the lightest Higgs boson. Unlike in MSSM, in the majority of 
models under study this state appears to have a significant component of the pseudoscalar 
A . We leave the details of the decay width for later, when more experimental information 
could become available. 
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FIG. 3. BRffff -> x°iXi) as a function of 8 S and tan/3, and BR(F 1 ° -> 66), BR(F 1 ° -> t+t - ), 
BR(if 1 ° WW*) and BR(i? 1 ° as functions of 9 S for the CP- violating versions of 

{7(1)77, [/(l)^, U(l)i, U(1)n and £7(1),/, models. When available, we also show the value of the 
corresponding SM quantity. 



V. DISCUSSION AND CONCLUSION 



The recent discovery of a Higgs-like boson at the LHC does not preclude the possibility 
of beyond the Standard Model (BSM) physics. With increased energy and luminosity, the 
couplings of the Higgs boson to SM particles will be measured with increased precision. 
In addition to the SM modes, the BSM Higgs boson can decay invisibly (to neutralinos, 
heavy neutrinos, or additional scalars). Our work investigates such a possibility, in a U(l)'- 
extended supersymmetric model, by analyzing the decay patterns of the lightest neutral 
Higgs boson. This study is motivated by the fact that the composition of the Higgs bosons 



25 




FIG. 4. M H o as a function of 9 S (the phase of the new singlet S) and 8f (the phase of the soft 
coupling At) for the CP-violating versions of the ?7(1)„, U(l)s, U(l)j, U(1)n an d U{\)ip models. 



is different from one in the SM or MSSM and hence, production and decay mechanisms are 
affected. Also significant is that U(l)' models, unlike the SM, predict a light Higgs boson 
(m H o ~ 125 GeV) naturally. 

We chose anomaly-free versions of U(l)' motivated by breaking of string-inspired E e SSM, 
and study the effects for both the CP-conserving and CP-violating scenarios, and compare 
the lightest Higgs boson production and decay to that in the SM. Our analysis has two goals: 
one is to analyze effects of CP violation on Higgs masses and decays, the other is look for 
differences among each of the U(l)' models for decay patterns, and identify characteristic 
signatures. 

We perform a complete study of Higgs sector of the effective U(l)' models, starting with 
calculation of masses and mixings in the Higgs sector, and including corrections from the 
stop and sbottom sector to one-loop level. Then we introduce benchmark scenarios for each 
EqSSM motivated U(l)' model, defined in terms of soft parameters, and the Higgs, Z' and 
sparticle spectra obtained for the benchmarks. We include a complete spectrum for the 
neutralinos, and include the saturation of the relic density constraint for each of the five 
versions of the U(l)' models. Our mass spectra calculation is restricted by the inclusion of all 
the known constraints on the low energy spectrum, and including all the recent constraints 
on the lightest Higgs boson mass, and also for rare decays and cosmological constraints. 

We then investigate the cross sections in channels (the vector fusion channel and the 
associated Higgs production with a vector boson) most propitious to look for the Higgs 
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boson to decay invisibly. While the cross sections are not significantly affected by the CP 
phases (coming from the effective \i parameter and the scalar trilinear couplings), the masses 
and the branching ratios show significant variations. With one exception, the decay into the 
lightest neutralino pair is significant in all, and dominant in two of the five U(l)' models 
under investigation. The invisible decay comes with, sometimes a suppression of the bb 
decay mode, from 60% to as low as 36%, except for U(l)s and U(l)i, where the branching 
ratio is enhanced with respect to the SM, up to 73%, in the absence of CP violating phases. 
All models exhibit a strong suppression of r + r~ mode (by a factor of 2-3), of WW* (by 
a factor of 2-5) and of 4£ by the same factor. Some of these branching ratios seem to 
be in agreement with the present LHC data although the measurements are not 

yet precise enough for a conclusive statement. The strong suppression in all U(l)' models 
of the decay into WW* can be traced to the mixing with the singlet, the pseudoscalar, 
and the CP-phase contribution, all which are known to modify the couplings of the Higgs 
boson with respect to their SM values. Overall, we find that Higgs phenomenology in 
U(l)' model is significantly affected by the CP phases, especially 9 S , and yields distinct 
signatures. The resulting signatures are unlike those of the NMSSM with CP violation, 
where the branching ratios of the lightest neutral Higgs boson are fairly independent of 
the values of CP phase 9 S |47fl. Some of the signals in U(l)' are typical of the anomaly- 
free versions of the models studied, others are characteristic for a scenario (such as the 
enhancement of the branching ratio into bb in U(l)s and C/(l)j). While other generic tests 
of CP violation in the supersymmetric sector exist, such as measuring chargino polarization 
{20 1 , the dependence of the masses and decay patterns of the Higgs boson with the phases 
are a much more promising indications for CP violation in U(l)'. Such signatures can be 
probed at the LHC, and are within reach at y/s = 14 TeV with luminosity C = 100 fb^ 1 . 

The decay patterns would enable to distinguish U(l)' models from the SM, but also from 
each other. For instance, U(l)s and U(l)i show some similar decay patterns, insofar as the 
decay Hf — > bb is dominant. Among all the models studied, U(l)s is the only one where the 
branching ratio of Higgs decay into neutralinos is below 10%; while in U(l)i the branching 
ratio into invisible modes is in the 10-20% range. In U(l) v , U(1)n and U(l)^, the partial 
width into the invisible mode is significant, but in U(l) v it is still slightly below that into bb. 
Distinguishing between U(1)n and U(l)^ could also be based on the branching ratio into 
the invisible channel, which can be over 50% in U(l)^, but under 50% in U(1)n. 
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The characteristic signatures at the LHC would be distinctive kinematic distribution of 
the two quark jets in the Higgs production through vector boson fusion, compared to the Zjj 
and Wjj backgrounds. In the Higgs production with an associated vector boson, the ZH 
associated production seems more promising, as a clean signal in the dilepton +$t channel 
will have little background, unlike the WH model where the single lepton +$t suffers from 
large background effects from off-shell Drell-Yan production, as previously discussed in the 



literature 



This scenario also has consequences for other neutral Higgs states, and for 



the charged Higgs, the analyses of which await more data. 
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Appendix A: Explicit Mass Formula 



In these appendices we give the complete and detailed analytical expressions used in our 
calculations. 



1. Scalar Top and Scalar Bottom Masses 



Stop and sbottom mass-squared matrices show clearly the differences between the MSSM 
and U(l)' extended models. As can be seen from the following expressions, extra charges 
and gauge couplings affect LL and RR entries especially if the vacuum expectation value of 
the S field is sizable (v$ > 1 TeV). 

The entries of the field dependent M 2 for scalar top are given by 

M\ L = M 2 ~ + Y 2 \H U \ 2 - - A {gl - 9 f){\H u \ 2 - \H d \ 2 ) + g 2 Y ,Q Q {Q u \H u \ 2 + Q d \H d \ 2 + Q S \S\ 2 ) , 

M 2 RR = M 2 D + Y 2 \H U \ 2 - ^f(\H u \ 2 - \H d \ 2 ) + g 2 Y ,Qu{Q u \H u \ 2 + Q d \H d \ 2 + Q S \S\ 2 ) , 

m lr = M 2 R [ = Y t (A*tH° u * - Y s SH° d ), (Al) 
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similarly for the scalar bottom mass-squared, we have 

M 2 LL = Af| + Y b 2 \H d \ 2 + \{gl + ^-)(\H U \ 2 - \H d \ 2 ) + g 2 Y ,Q Q {Q u \H u \ 2 + Q d \H d \ 2 + Q S \S\ 2 ) , 

M 2 RR = M 2 ~ + Y 2 \H U \ 2 + ^{\H U \ 2 - \H d \ 2 ) + g 2 Y ,Q D (Q u \H u \ 2 + Q d \H d \ 2 + Q S \S\ 2 ) , 

m lr = M 2 R [ = Y b (A* b H° d * - Y S SH° U ). (A2) 



2. Neutral Higgs Boson Masses 



The neutral Higgs masses are obtained by diagonalizing the 4x4 matrix in Eq. ([TBI) . 
The explicit values of the entries are: 



M 2 



11 



2 2 
77i~ /77, ~ 

S,(3A?(A2( 2 y>3 l n (_^) + ^^(A 6 C 6 F 6 y 6 2 + A t C t F t Y 2 )) 

ml 



+ 2/i 2 G b r f X(AC^ - fiVuY) + 6A 2 t A 2 b G t Y t 4 v u (fiC t v d - A t v u ) 2 + 647r 2 A 2 A 2 A^ 



+ 12A t A 2 (G t - 2)Y t i Ay u (^C t v d - A t v u ) + fi X ^Pd^t 



(A3) 



M\ 2 



3S^,S t 



A2s i (E 6 (3/iA b F fe 2 (2G b F f) 2 (A b C b ^ - fiv u )(A b v d - ftC b v u ) + C b A 2 b F b 



- 32vr 2 A 2 t Wl A ud ) + 6^A 2 (G b - 2)Y b 4 v d (/iv u - A b C b v d ) + fi X A 2 b E b ) 
+ 6vA 2 b Z b Y t 4 {fiv d - A t C t v u ){A t G t E t {fxC t v d - A t v u ) + (G t - 2)A t 2 <) 



+ 3fiA t A 2 b X b C t F t Y t 2 A 2 X t 



(A4) 



.M 2 



13 



3f S S t 



A 2 S t (A 2 (3/iF fe y fe 2 (2^ u - A 6 C6V d ) + 32tt 2 ^ u Aus) 



+ 6/2 2 G b Y b % u (A b C b v d - /jv u y - nx&iv d ) - 3fiA t A 2 b C t v d F t Y t 2 A 2 t E t 
- 6fiA 2 b v d Y t \fiv d - A t C t v u )(A t G t Z t (iMC t v d - A t v u ) + (G t - 2)A 2 v u ) 



(A5) 



M 2 U 



Z t (fiA b G b Y b 4 S b A 2 t ^v u - A b C b v d ) + A 2 t A 2 b G t Y?S t (A t v u - fxC t v d )) 



A t A 2 b {G t -2)Y t A S t A 2 v u 



(A6) 
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M 2 



K 
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3T, b v d 



1 2 

mf m~ 

j: b (3A 2 (A 2 b (2Y b 4 v d H^^) + ^v u (A b C b F b Y 2 + A t C t F t Y 2 )) 

Tfly, 



+ 2A 2 b G b Y b % d (A b v d - fiC b v u ) 2 ) + 6pi 2 Aiv d G t Y t *(piv d - A t C t v u ) 2 + 64vr^A^^A, 2 ) 



l2A b k 2 b (G b - 2)Y b % 2 A 2 t (A b v d - ^C b v u ) + // X A^£ 6 A t V 



(A7) 



K 



Mi, =- 

23 3X b v s 



A 2 t {Y, b {3 l iA 2 {F t Y?{2 l iv d - A t C t v u ) - A b C b F b Y 2 v u ) 

- 6fiA b G b Y b v u (fiv u - A b C b v d )(fiC b v u - A b v d ) + 327r 2 A^ (i A ds t>§) 

+ 6/iA 2 (G 6 - 2)Y b 4 v d v u (A b C b v d - ptv u )) + 6pL 2 A 2 b Z b v d G t Y t 4 (piv d - A t C t v u ) 2 



(A8) 



, y9 2kuuj 
M 2 24 = —^ 

^bV S 



E b (AiG b Y 4 S b Af(A b v d - ptC b v u ) + fiA t A b G t Y 4 S t (fj,v d - A t C t v u )) 



A b A 2 b {G b -2)Y b 4 S b v d A\ 



(A9) 



M 



33 " 3v? 



A 2 (3^v u (A 2 b v d (A b C b F b Y 2 + A t C t F t Y 2 ) + 2^G b Y b 4 v u (A b C b v d - fiv u f 



+ Q4ir 2 A 2 b \ s v 4 ) + Q^A 2 b v 2 d G t Y 4 ^v d - A t C t v u f + // X A^A t X), 



(A10) 
(All) 



M 



2K,[1 2 U! 



34 



A b G b Y 4 S b A 2 t v u (piv u - A b C b v d ) + A t A 2 b v d G t Y 4 S t (pLv d - A t C t v u ) 



, (A12) 



M 2 M = 



3v d v 2 s v u 



3A 2 t (A 2 b (A b C b F b Y b 2 + A t C t F t Y 2 ) + 2 l iA 2 b G b Y 4 S 2 b v d v u 



+ ^AlAlvdGtYfS'tvu + X A^A 



4 o2„ 



2 A 2 



(A13) 



3. CP-odd Tadpole Terms 



Explicit form of the CP-odd tadpole terms are 

T 4 = MsUdSin(fe + S ) + ^^v d (A b F b Y 2 S b + A t F t Y 2 S t ), 
% = fiA s v u sin(0 s + 9 S ) + -^3fxv u (A b F b Y 2 S b + A t F t Y 2 S t ), 



32tt 2 

fj,A s v d v u sin(0 E + 6>s) 1 

/6 — T 



v 5 



32tt 2 vs 



3fiv d v u (A b F b Y 2 S b + A t F t Y 2 S t ). 



(A14) 
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4. Charged Higgs Boson Masses 



Finally, the charged Higgs mass is obtained by diagonalizing the matrix in Eq. (l2~Uj) . One 
of the eigenvalues will be the Goldstone boson needed to give mass to the W ± boson, the 
other is the real charged Higgs mass. The explicit entries in (120 \ are: 



M 



2± 
11 



3v 2 viX t v u 



K Aiv d A 2 t (E t (fivy s (3A b C b F b Y b 2 + x) + 3fiA b F b Y b 2 S b v 2 s n 



2„,2 
u 



+ v d v u (87r 2 v 2 d ( g y s - V) - 3^ 2 F b Y 2 v 2 s )) + 3Y t 2 v 2 (A t F t J: t (v u (fiS t v u - A t v d ) + fiC t vj) 



mi 



v d Ztv u (F t + G t -2) + v d (G t - 2)A 2 t v u ) + 6v d Y t ^Z t vZ(ln(-^) - 1)) 



(A15) 



M 2± - - 

Mw ~ 3v 2 E b v 2 



K A 2 b A 2 t {E b {3fxA t F t Y t 2 v 2 s {C t vl + v 2 d S t ) + v u {-3^ 2 v d F t Y 2 v\ 



+ 8?r v d v u (g 2 v s - 4/i ) + ^x v s v u)) + 3F b v s (A b F b E b (v u (iiC b v u - A b v d ) + ^iS b v d ) 



mi 



- Z 2 b v d v u (F b + G b -2) + A 2 (G b - 2)v d v u ) + 6F b 4 S^>^ n (ln(^f ) - 1)) 



(A16) 



M 



2± 
21 



1 



3v 2 vlE t 



KAiA 2 (Z t (ij,vy s (3A b C b F b Y 2 + x ) + 3fiA b F b Y 2 S b v 2 s v 2 



+ v d v u (87r 2 v 2 d (g 2 v 2 s - A/2 2 ) - 3fi 2 F b Y 2 v 2 s )) + 3Y 2 v 2 s (A t F t E t (v u (fiS t v u - A t v d ) + fiC t v 2 d ) 
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v d Y? t v u {F t + G t -2)+ v d (G t - 2)A 2 v u ) + 6^F>|S t ^(ln(^|) - 1)) 



(A17) 



1 



M ™ ~ 3v 2 E b v d v 2 



KA 2 b A 2 t v u (i: b (3fiA t F t Y t 2 v 2 (C t v 2 u + v 2 d S t ) + v u (-3fi 2 v d F t Y 2 v 2 s 
+ 8n 2 v d vl(gjv 2 s - 4/i 2 ) + nxv 2 s v u )) + 3Y b 2 v 2 s (A b F b T lb (v u (fiC b v u - A b v d ) + fxS b vj) 



- E 2 b v d v u (F b + G b -2) + A b (G b - 2)v d v u ) + 6Y*Z b vy s v u (\n 



mi 



Q 2 



i)) 



(A18) 



5. Auxiliary Expressions 



In the above expressions, we use the following short-hand notations: 



X = J1024ttM| - 9(A b F b Y 2 S b + A t F t Y 2 S t ) 2 , 



(A19) 
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and 



A u d = QdQudY' ~ ~j~ + ^5 5 

Ads = QdQ5^ + ^ 2 ; 
A us = QsQ u ^+^ 2 . 

(A20) 
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